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Interacting electrons in quantum dots with large Thouless number g in the three classical random 
matrix symmetry classes are well-understood. When a specific type of spin-orbit coupling known 
to be dominant in two dimensional semiconductor quantum dots is introduced, we show that a new 
interacting quantum critical crossover energy scale emerges and low-energy quasiparticles generically 
have a decay width proportional to their energy. The low-energy physics of this system is an example 
of a universal interacting crossover regime. 

PACS numbers: 73.50.Jt 



The statistics of the single-particle states of meso- 
scopic systems with disorder or chaotic boundary 
scattering [1 J2| a re controlled by Random Matrix The- 
ory (RMT) 3], as long as the states are separated by less 
than the Thouless energy Et (related to the ergodiciza- 
tion time for a particle r erg by the Uncertainty Principle 
Et = h/T erg ). For mean single-particle level spacing 5, 
the Thouless number is g = Et/S. 

Since disorder breaks all the spatial symmetries, only 
time-reversal T and possibly Kramers degeneracy re- 
main. There are three classical symmetry classesQ, 
the gaussian orthogonal ensemble or GOE (T intact, 
no spin-orbit coupling), the unitary or GUE (T bro- 
ken), and the symplectic or GSE (T intact, with spin- 
orbit coupling). More recently, other classes have been 
identified for disordered superconductors Q and quan- 
tum dots constructed from two-dimensional semiconduc- 
tor heterostructures with spin-orbit coupling^. We will 
focus on a symmetry class in the latter case (which we call 
the Aleiner-Falko (AF) class Q) in which, after a canon- 
ical transformation, the single-particle S z is conserved, 
while S 2 is not||. 

The idea of a crossover between two symmetry classes 
will play a central role in this paper. Consider a system 
in which the Hamiltonian is crossing over from the GOE 
to the GUE0], acheived, e.g., by turning on the orbital 
effects of a magnetic field. For Thouless number g 1 
of the original GOE, the g x g crossover Hamiltonian is 



Hx(ot) — Hqoe H — ^Hgue 
V9 



(1) 



where a is the crossover parameter. Properties of 
eigenvector^, 0| correlations have been computed in the 
crossover. For g ^S> 2a 2 = gx ^> 1, the following 
ensemble-averaged correlations hold for the eigenstates 
ipfi(i), where /i^y label the states and k, I the orig- 
inal orthogonal labels: 
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The last term on the second line shows the extra corre- 
lations induced in the crossover 0- The crossover scale 
Ex = gxS/TT represents a window within which GUE 



correlations have spread, while GOE correlations remain 
at high energies. 

The crossover from the spin-rotation invariant GOE 
to the AF class is a GOE — > GUE crossover, where the 
"magnetic flux" has opposite signs for opposite eigenval- 
ues of S z || ■ If the linear size of the system is L and the 
spin-orbit scattering length is £ 3> L, this new AF sym- 
metry class manifests itself below Ex — (j) 4 Et- The 
crossover to the fully symplectic GSE occurs @ at the 
parametrically smaller energy scale of (j) 6 Et, set here 
to 0. 

Turning from single-particle physics to interactions, 
for small to moderate r s the "Universal Hamiltonian" [1, 
Q is known to contain all the relevant couplin gs [loll at 
low energies in the renormalization group sense [llj. 



Uo^ 2 
2 



JS 2 + XT^T (3) 



Here TV is the total particle number, S is the total 
spin, and T = X) c /94 c /9-T- has a charging energy 

U, an exchange energy J and a superconducting cou- 
pling A. This last term is absent in the GUE, while 
the exchange term disappears in the GSE. In the large- 
g limit only interaction terms which are invariant under 
the symmetries of the one-body Hamiltonian appear in 
0- At larger r s , the system enters a quantum 
critical regime[l2j connected with the impending Pom- 
ranchuk transition [3, [l4[ ■ 

In the GOE, ignoring the Cooper coupling (two- 
dimensional semiconductor quantum dots do not super- 
conduct) and tuning J one sees the mesoscopic Stoner 
effect @, H. Take for illustration an evenly spaced set of 
levels with level spacing 8. Since both S 2 and S z are con- 
served, we can focus on the ground state with S z = S, 
and find its energy for an even number of particles to 
be E gs (S) = S 2 S - JS(S + 1). Defining J = J/5 and 
minimizing with respect to S leads to steps from S = 
to S = 1 at J = |, from 5=1 to S = 2 at J = f 
etc. Including the mesoscopic (sample-to-sample) fluctu- 
ations of the energies and the matrix elements leads to 
probability distributions for spin 50, HI [l^]. Note that 
the S z = S ground state is electronically uncorrelated (a 
single Slater determinant, not a superposition). 
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FIG. 1: Phase diagram in the Ex — J plane showing the lines 
of ground state S z transitions. Note the smooth approach to 
the spin-rotation-invariant result as Ex — * 0. 



What happens when spin-orbit coupling of the AF 
class is introduced into a slightly generalized form of 



H 



y] e i-is c }is c fj.s 



J Z S Z j{Sx + Sy 



(4) 



Here the basis \i labels the eigenbasis of the single-particle 
AF crossover hamiltonian[5j , and s is the eigenvalue of 
S z = § Yl c It c mT — c l i c a4 • Sx,y are expressed fl6| in terms 



of the combinations S± = S x ± iS y as 



S 4 
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with S- = (S+y. The matrix M a p depends on the par- 
ticular realization, and has the ensemble average (l6) 

ExS/n 



{\M at 



Ex 2 + (m - £/3l) 2 



(6) 



At strong spin-orbit coupling!!, 0, [l6j], noting that 
-J{S 2 X + Si) is irrelevant in the RG sense [13, HE fl7j ] ■ 
we end up with the Hamiltonian of Eq. ^ with J = 0. 
The ground and excited states of this Hamiltonian have 
definite S z and are electronically uncorrelatcd. 

We are now ready to state our central results. In the 
crossover to the AF class, a quantum critical regime[l2j 
emerges at a many-body quantum critical crossover scale 
Eqcx = fl — J) Ex (a result the author obtained 
previously [17[ in the restricted case J z =0). In contrast 
to the limits of zero and strong spin-orbit coupling, in the 
crossover the ground and low-lying states are electron- 
ically strongly correlated. Transverse spin fluctuations 
have a nonzero density at low energies and fermionic 
quasiparticles become very broad at low energies (as long 



as ground state S z ^ 0). Finally, the mesoscopic Stoner 
effect is smoothly pushed to higher J as spin-orbit cou- 
pling increases. 

We will set J z = J in Eq. (|4|) henceforth, since that is 
the correct starting point for a Hamiltonian deep in the 
Thouless bandflol. Iig. [lij . We decompose both the inter- 
action terms by introducing the Hubbard-Stratanovich 
fields h(t), q{t) = X(t) + iY(t), and q*(t) to get the 
T = imaginary time action 



A= J dt ^as[(dt + e a )S ss > - ^(o- z ) ss >]ip as > 

— og a>,ss' 
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X(t),Y(t) are fluctuating fields and are integrated out, 
but h(t) acquires an expectation value (also called h). For 
our picket fence spectrum e„ = 5(n — j) with chemical 
potential at 0, at h = all states n < are occupied 
while all states n > 1 are empty When h lies between 
h n = 5(2n — 1) and h n+ i = <5(2n+l), at T = the single- 
particle states between — n + I and n singly occupied, 
and S z = n = S. We will find saddle points for h in each 
of the intervals h„ < h < h n+ \, thereby obtaining the 
ground state energy as a function of S = n, and look for 
the lowest one. 

First integrate out the fermion fields and obtain a 
quadratic effective action for h, X, and Y. 

Ae ff = jiy+ wi y M)|, (i-w) 

— oo 

_| 4j — L Jxi(iu) (8) 

The cross terms are a consequence of [S x , S y ] = iS z , and 
Xr and \i are t ne rea l an d imaginary parts of x> the 
fermionic transverse spin-susceptibility. 



}^ \M mn \ 



> N F (e m ^) - N F (e nl ) 
iuj + e n ; - e r „ T 
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where Np is the Fermi occupation. To make further 
progress, we assume that Ex S, which allows us to 
convert the sums over states into integrals, and also re- 
place the sample-specific value of \M mn \ 2 by its ensem- 
ble average. Such self-averaging occurs naturally in the 
large- N limit (H, [lij . The dominant contribution to % is 



1 E x - iE s sgn(o;) 
(5 |w| + Ex — ih sgn(tu) 



(10) 



where E$ = 2SS and sgn(uj) is the sign of u>. Next, 
the integration over X and Y results in a fluctuation 
contribution to the effective action 
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where the argument of the logarithm is the determinant 
of the matrix of the quadratic form of X and Y. Afi uc 
is logarithmically divergent due to a limitation of the 
ensemble averages Eqs. (|2I6[) for large energy separations 



Ex > Et- Cutting it off, discarding terms of the form 
\og(E T /E x ) which are independent of S, defining Ah = 
h — JEs , and adding the one-body energy we find 



Ae ff (S,h) = S 2 6 + ^--hS+- 

4 J 7T 



Ah tan 



Ah 
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Eq. (JT2J) is one of the central results of this paper. Note 
that the term in the brackets is or order 1/5 compared 
to the first three terms (h will turn out to be order S). 
From it we find that the saddle point value of ho satisfies 



fraction of S and will be ignored below. Going over to 
the retarded commutator with the standard replacement 
ioj — ► lj + i0 + , we find the spectral function for trans- 
verse spin excitations 



~ 2J , _ x Ah 
ho = JE S H ( tan — 

7T EqcX 



tan" 1 |L) (13) 



We must also ensure that h n < ho < h n +\ for S = n. 
Using h = ho in Eq. Q12[) gives the ground state energy 
for that S. Even though this result has been derived 
for Ex 3> S and S 3> 1, note that one recovers the 
correct spin-rotation-invariant result for all S on taking 
the Ex — * limit first in Eq. (|12p and solving it to 
obtain the ground state energy. Fig. [T] shows regions 
in the Ex, J plane with different ground state S z = S. 
Note the smooth approach to the spin-rotation invariant 
results as Ex — ► 0. 

Now consider the low-lying excitations, which are 
bosonic spin excitations and fermionic quasiparticles. 
The transverse spin correlator 

V(t) = -<T t 5+(t)5_(0)> = -±.(T t <f(t)q(0)) + - 

(14) 



function of Matsubara frequency (at T = 0) is 



V(iuj) 
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Eqcx + iAh sgn(u) 
The last term in the denominator is never more than a 



B(u) = -2Im(D ret (uj)) 
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For S = this reproduces the scaling function com- 
puted earlier [It]]. More importantly, for 5^0 there is 



a nonzero density of spin excitations even as 



< 1 



(but u>if). These excitations are related to the 2S+ 1- 
degenerate ground states of the spin-rotation-invariant 
system (since S 2 is not conserved here, the ground state 
contains a superposition of many values of S 2 ). These 
low-energy excitations make the system strongly corre- 
lated in the electronic sense, that is, the ground and low- 
lying states are no longer described even approximately 
by single Slater determinants. 

The decay of an electron of energy e to leading order 
occurs via the emission of a single spin excitation. Using 
the Fermi Golden Rule, we get the decay rate 



r(e) ~ 2ttJ 2 / de'\M(e - e')\ 2 B{e - s')p(e') (17) 



where p{e') is the electronic density of states at energy e' 
and M is the matrix of Eq. ©. Using p(e') = 1/8, and 
Eqs. (|6I16[) we obtain 
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This is the other central result of this paper. At low en- 
ergies e <C Eqcx, P, \M\ 2 , and B are constant, leading 
to a decay rate which goes as T{e) ~ Ex^Eqcx £ ' w ^ c ^ 
can exceed e for J — * 1, leading to ill-defined quasipar- 
ticles. At high energies e 3> Ex (but e -C Et), the 



decay rate due to spin excitations goes to a constant. 
There will be an additional Fermi liquid broadening fl9| . 
r ~ e 2 /Et due to neglected interactions. Thus we have 
an unusual situation in which the quasiparticles are bet- 
ter defined at high energies than at low energies, because 
the high-energy physics is controlled by the weakly inter- 



4 



acting spin- rotation- invariant Hjj 0, Q . 

In summary, we have established that a crossover from 
the spin-rotation-invariant GOE to the AF class Q in 
two-dimensional semiconductor quantum dots makes the 
ground and low-lying states of the system strongly cor- 
related in the electronic sense, with a nonzero density of 
states for transverse spin excitations at very low ener- 
gies, and a decay rate T(e) ~ e for low-energy fermionic 
quasiparticles (for ground state S z ^ 0). The results 
are universal since all energies are <C Ex, and apply to 
disordered as well as ballistic/chaotic dots. 

Random matrix crossovers thus offer us access to uni- 
versal interacting crossover reaimesfv^ . Such regimes 
are dominated by many-body correlations and are dis- 
tinct from single-particle RMT ensembles 0, 0, and 
offer the possibility of tuning the many-body quantum 
crossover scale Eqcx by varying the single-particle Ex- 

Another example of such a universal interacting 
crossover regime is a superconducting dot verti cally 
tunnel-coupled to a normal dot with an orbital flux[2fj. 
The crossover scale is varied by varying the tunneling 
strength. In our problem one can imagine a pair of verti- 
cally coupled quantum dots with one of them fabricated 
of a different material having a much larger spin-orbit 
coupling than the other. Ex (and thus Eqcx) would be 
tuned by varying the tunneling strength. 

Let us turn to some caveats. We have focused on an 
even number of electrons in the dot but expect simi- 



lar results for an odd number. Our large-iV approach 
is valid for Ex ^> S and 5 » 1. However, even in this 
regime, in a tiny region of width A J ~ <5 2 / (ExS) around 
the transition between S z steps, the value of ho is such 
that the fermionic susceptibility exceeds 1/J, and a small 
saddle-point value for |g| will be generated, leading to a 
more complicated effective action. Apart from poten- 
tially missing effects in this region, the large- N approach 
seems to capture most of the physics. We have assumed 
a spectrum with equal level spacing, and the ensemble 
averaged form for the matrix elements M a p. In real dots 
in the crossover, both fluctuate[||, @, One thus ex- 
pects mesoscopic (sample-to-sample) fluctuations in the 
ground state value of S z = 50, M, EE EH leading to a 
probability distribution for S given J and Ex ■ Our meth- 
ods can be used to complement exact diagonalizations pi) 
in this case. 

Many open questions remain, such as how to charac- 
terize the state at and near the transition between S z 
steps, how the state responds to an in-plane B field, 
whether the strong electronic correlation has any signa- 
tures in the zero-bias conductance, and whether one can 
classify universal interacting crossover regimes into uni- 
versality classes. The author hopes to explore these and 
other questions in future work. 

It is a pleasure to thank R. Shankar and and especially 
Yoram Alhassid for comments on the manusript. 
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